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We point out that domain wall formation is a more common phenomenon in the Axiverse than previously 
thought. Level crossing could take place if there is a mixing between axions, and if some of the axions acquire 
a non-zero mass through non-perturbative effects as the corresponding gauge interactions become strong. The 
axion potential changes significantly during the level crossing, which affects the axion dynamics in various ways. 
We find that, if there is a mild hierarchy in the decay constants, the axion starts to run along the valley of the 
potential, passing through many crests and troughs, until it gets trapped in one of the minima; the axion roulette. 
The axion dynamics exhibits a chaotic behavior during the oscillations, and which minimum the axion is finally 
stabilized is highly sensitive to the initial misalignment angle. Therefore, the axion roulette is considered to be 
accompanied by domain wall formation. The cosmological domain wall problem can be avoided by introducing 
a small bias between the vacua. We discuss cosmological implications of the domain wall annihilation for 
baryogenesis and future gravitational wave experiments. 


Introduction - There may be many axions in nature. In 
string theory, compactihcations of the extra dimensions often 
offer a large number of axions [1, 2]. The Universe with a 
plenitude of such axions whose masses range over many or¬ 
ders of magnitude is called the Axiverse [3]. The striking fea¬ 
ture of the axions is that they respect shift symmetry, which is 
unbroken at the perturbative level. Non-perturbative effects, 
however, typically break the shift symmetry to a discrete one, 
generating a periodic potential for the axions. Some of the 
axions may remain relatively light and play an important role 
in cosmology such as inflation [4, 5], dark energy [3, 6], dark 
matter [7, 8], and baryogenesis [9-12]. 

The axion mass is not necessarily a constant, but can vary 
with time. For example, the QCD axion acquires a non-zero 
mass at the QCD phase transition, while it is massless at high 
temperature. Similarly, some of the axions may acquire a 
mass after inflation. This is expected to be the case if the the 
hidden sector is heated to high temperatures where the hidden 
gauge coupling is weak. As the temperature of the hidden sec¬ 
tor decreases in the course of the cosmic expansion, the axion 
mass is generated when the hidden gauge interactions become 
non-perturbative. 

In general, axions have mass and kinetic mixings with each 
other [13]. The level crossing takes place if the mass of one 
(or more) of the axions increases with time due to the non- 
perturbative effects and becomes heavier than other axions. 
During the level crossing, the axion potential changes signifi¬ 
cantly, which affects the axion dynamics in various ways. The 
level crossing phenomenon among the axions in the Axiverse 
has not attracted attention so far to the best of our knowledge, 
and we shall study its cosmological impact in this Letter. 

The axion starts to oscillate when its temperature- 
dependent mass becomes comparable to the Hubble param¬ 
eter. If the axion starts to oscillate well before the level cross¬ 
ing, the change of the axion potential is sufficiently slow com¬ 
pared to the oscillation period. In this case, we find that the 
resonant transition takes place at the level crossing a la the 
MSW effect in neutrino physics [15, 16]. On the other hand. 


if the commencement of the axion oscillations is close to the 
level crossing, the axion does not really oscillate along one 
direction, but starts to run along the valley of the potential, 
passing through many crests and troughs, until it gets trapped 
in one of the potential minima. We find that which minimum 
the axion is finally stabilized is highly sensitive to the initial 
misalignment angle, and as we shall see, its dependence is 
rather chaotic. We call the chaotic run-away behavior the “ax¬ 
ion roulette.” 

Because of the high sensitivity to the initial position and 
possible spatial instabilities, the axion roulette is expected to 
be accompanied by domain wall formation. Therefore, do¬ 
main wall formation might be a more common phenomenon 
in the Axiverse than previously thought. We will determine 
the conditions for the axion roulette to take place. The domain 
walls are cosmologically problematic, and they must be either 
inflated away, or unstable and decay rapidly. We shall discuss 
the implications of the domain wall collapse for baryogenesis 
and future gravitational wave experiments. 

Lastly let us briefly mention the works in the past. The res¬ 
onant transition between the QCD axion and the hidden axion 
was studied in the pioneering paper by Hill and Ross [17]. 
Recently two of the present authors (NK and FT) studied the 
resonant transition in Ref. [18], where the adiabatic invariant 
was correctly identified, and it was shown that not only the 
axion abundance but also its isocurvature perturbations can be 
significantly suppressed if the adiabaticity is weakly broken 
by the initial condition of the axion close to the hilltop. As 
emphasized in Ref. [18], the level crossing phenomenon is not 
limited to the QCD axion and the hidden axion, but it occurs 
commonly between any axions with a mixing. The purpose of 
this Letter is to show that the level crossing leads to not only 
resonant transition, but also domain wall formation. 

Level crossing - Let us consider the low-energy effective 
Lagrangian for axions ai and a^', 

£ = ^ - 14(01,02) - 1^2(01,02) ( 1 ) 


2 


with 


V"i(01,02) = Ai ( 1 - cos ( + 712 ^ 


^2(01,02) = maiTYfi ( 1 - cos ( ^ 


where ma{T) is given by 
ma{T) = min 


T 

Wa 1 ^ ) , Wfl 


( 2 ) 

( 3 ) 

( 4 ) 


Here we have defined rua = AI//2 and ni (712) and /i (/2) 
are, respectively, the domain wall number and the decay con¬ 
stant for oi (02). p is a negative constant, which depends on 
the details of the interactions responsible for V2. In the follow¬ 
ing numerical calculations we assume the radiation dominated 
Universe and set p = —3 for simplicity, but its precise value 
is not relevant for our main results. Note that the two axions 
have a mixing through Vi, and so, ma{T) is not necessarily 
equal to one of the mass eigenvalues. 

Let us define two effective decay constants for later use. 
When the curvature of Ui is larger than that of U2, one can 
define the effective decay constants / and F along the val¬ 
ley of Ui and the direction orthogonal to it, respectively. 
Specifically, the lighter direction is determined by niaij fx + 
^2027/2 = const.. The effective decay constants are given 
by 
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and one can see that, for large rii and/or 712, there is a hier¬ 
archy between / and F. Such hierarchical decay constant at¬ 
tracted attention in a context of the so called alignment mech¬ 
anism for large-field inflation [ 19 - 21 ]. The reason for this 
choice of the axion potential will be clear shortly. 

We assume that, while the potential height of Vi is constant 
with time, the potential height of V2 grows with time. This is 
the case e.g., if the gauge sector responsible for V2 is heated 
to high temperature where the gauge interactions are weak. 
Then, at high temperature, there is a massless flat direction 
along 7 iiai//i + 772027/2 = const.. As the temperature de¬ 
creases, ma{T) turns on and the height of V2 increases. We 
assume that the mass of 02 finally becomes heavier than oi, 
and then it becomes constant in the low energy. At low tem¬ 
peratures, the heavier mass eigenstate is almost 02, while the 
lighter one is Oi. See Fig. 1 where we show typical evolu¬ 
tion of the two mass eigenvalues, (777^,777//) with rriL < rnn, 
while the height of V2 increases due to non-perturbative ef¬ 
fects. Here the mass eigenvalues are evaluated at the potential 
minimum. Here and in what follows we adopt 771 =772 = 10 , 
/i = /2 = 10 ^® GeV and {K1/K2Y = 0.02 as reference 
values. The level crossing occurs at ruat ~ 10 when the two 
mass eigenvalues become comparable. 



FIG. 1. The time evolution of the mass eigenvalues of two axions, 
mu and 777 _l. 


The axion potential changes significantly during the level 
crossing. This might affect the axion dynamics if the ax¬ 
ion has already started to oscillate before or around the level 
crossing. If the axion starts oscillations along the valley of Ui 
much before the level crossing, the change of the axion po¬ 
tential is so slow compared to the typical oscillation period 
that the axion number density in the comoving volume be¬ 
comes the adiabatic invariant. Then, resonant transition takes 
place during the level crossing a la the MSW effect [ 15 , 16 ]. 
As a result, the final axion density is suppressed by the mass 
ratio, mj^/rriH, compared to the case without the resonant 
transition [ 18 ]. The adiabaticity is violated if the oscillation 
amplitude is so large that the anharmonic effect is relevant. 
In this case, the resonant transition is not complete and both 
heavy and light axions are generated. Then, isocurvature per¬ 
turbations can be suppressed for certain parameters, as the 
produced heavy axions partially cancel the original fluctua¬ 
tions [ 18 ]. 

Axion roulette and domain wall formation - Now let us 
consider the case in which the axion starts to oscillate around 
or slightly before the time of the level crossing. Then, the 
scalar potential significantly changes even during the first few 
oscillations of the axion. As a result, the direction of oscil¬ 
lations is changed each time the axion oscillates. This is be¬ 
cause the axion effectively gets kicked into different directions 
around the turning points of oscillations, and this lasts until 
the level crossing is completed. Note that there is no adiabatic 
invariant in this process because the motion is not periodic at 
all. 

We have assumed that V2 eventually dominates over the 
scalar potential. Namely, at low temperatures T A2, the 
heavier mass eigenstate is approximately given by 02, while 
the lighter one is oi. Therefore, after the commencement of 
oscillations, the axion dynamics is likely confined in one of 
the valleys of V2, e.g., around 02 « 0 . On the other hand, 
it is not certain if the axion is quickly trapped in one of the 
minima of Vi. In fact, if the initial axion energy is sufficiently 
large, the axion will start to run along one of the valleys of 
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FIG. 2. The trajectory of two axions in the ai-a 2 plane. The ax- 
ion evolves from left to right. The region closed hy dashed-green 
(solid-orange) contour contains the potential minimum (maximum) 
atm4Tffi=Ai 

V2 soon after the onset of oscillation, passing through crests 
and troughs of Vi, until it gets trapped in one of the potential 
minima. 

We have numerically confirmed such a behavior. In Fig. 2 
we show the trajectory of the two axions in the 01-02 plane. 
We adopt the initial condition a2,i/fi — 1-5 and n\a\^il -F 
n2a2,i/f2 — 0 which minimizes the Vi, but this does not af¬ 
fect our main results. One can see that the axion first starts to 
evolve along the flat direction of Vi, and then, gradually starts 
to run toward positive values of oi along the valley of V2 (i.e., 
02 = 0 ). In this example, the axion is stabilized after passing 
through ten crests of Vi shown by the solid-orange lines. We 
call the run-away behavior of the axion dynamics the “axion 
roulette!’ 

The axion roulette takes place even if the axion mass is var¬ 
ied by about one order of magnitude, as long as the axion 
starts to oscillate slightly before or around the level crossing. 
According to our numerical calculations, the axion roulette 
works well for H\c/Hosc = 0 ( 0.1 — 1), where the subscripts 
Tc’ and ’osc’ imply that the variable is evaluated at the level 
crossing and at the commencement of oscillations, respec¬ 
tively. On the other hand, if the axion starts to oscillate well 
after the level crossing, the axion initially starts to oscillate 
along 02, and no motion along oi is induced, suppressing the 
chaotic run-away behavior. 

In order to understand the condition for the axion roulette 
to take place, let us estimate the energy density of the axion 
Pa,osc at the commencement of oscillations. In order to climb 
over the potential barrier of Vi, the initial oscillation energy 
must be greater than the height of Vi: 

Pa,osc ~ tJl^L.oscf^ ^ ~ ( 7 ) 

In deriving the above condition ( 7 ) we have assumed that the 
axion starts to move along the valley of Vi, because the cur¬ 
vature is still dominated by Vi until the level crossing, even if 
the potential height of V2 is larger than Vi. One can imagine 
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FIG. 3. The final value of (riiOi / fi +n 2 a 2 //2)/27r as a function of 
the initial misalignment angle of a 2 . The bottom panel is the zoom- 
up of the top panel between 02 ,i //2 = 1.5 - 1.5 -F lO”'^. 


a cosine potential V2 with grooves induced by Vi. Around the 
level crossing, the mass eigenvalues at the potential minima 
are comparable to each other. Therefore, if the effective de¬ 
cay constant along the flat direction is enhanced,i.e., f ^ F, 
the initial axion energy easily exceeds the height of the poten¬ 
tial barrier. This is the reason why we have chosen the axion 
potential ( 2 ) and ( 3 ), where the alignment mechanism is op¬ 
erative. Specifically, if one sets /i = /2 and ni = n2 = N, 
Eq. ( 7 ) yields N > ttih/itilIosc- Equivalently, one may in¬ 
troduce a hierarchy in the two decay constants from the be¬ 
ginning. In contrast to the application to large-held inflation, 
we do not need very large hierarchy in the decay constants or 
super-Planckian ones. The hierarchy by a factor of several to 
ten is sufficient. 

So far we have studied the spatially homogeneous evolution 
of the axion, neglecting spatial inhomogeneities. The axions 
acquire quantum fluctuations of order i 7 inf/ 27 r during infla¬ 
tion, if they are lighter than the Hubble parameter, iTinf ■ If the 
final minimum where the axion is stabilized is highly sensitive 
to the initial position, the domain walls are likely produced. In 
Eig. 3 we show the final value of (riiai/fi + n 2 a 2 // 2 )/ 27 r 
as a function of the initial misalignment angle of 02. The 
chaotic behavior becomes strong at large misalignment angles 
because the initial oscillation energy is increased. For larger 
N, the chaotic behavior is observed at smaller misalignment 
angles. Therefore the axion can be stabilized in different po¬ 
tential minima in each Hubble patch, even if the quantum fluc¬ 
tuation is much smaller than the decay constant. In addition, 
even at subhorizon scales, spatial inhomogeneities are likely 
to grow during the time when the axion runs through crests 
and troughs of Vi. Therefore, the axion roulette is most prob¬ 
ably accompanied by domain wall formation. 

Discussion - The domain walls are cosmologically prob¬ 
lematic, and they must be either inflated away, or unstable 
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and decay quickly before dominating the Universe. The latter 
can be realized if one introduces another small shift symmetry 
breaking term to lift the degeneracy of the potential minima. 
The domain wall collapse and the emitted gravitational waves 
have been studied extensively in the literature. [22-25] Our 
analysis can be straightforwardly applied to the mixing be¬ 
tween the QCD and hidden axions. In this case, the peak fre¬ 
quency of the gravitational waves is comparable to or lower 
than 0(10“^) Hz, which may be within the sensitivity reach 
of the pulsar timing experiments. [26, 27] 

So far we have introduced a large domain wall number 
ni = n 2 = N > 1 to realize the required hierarchy in the 
effective decay constants. Instead, we may introduce hierar¬ 
chy in the decay constants from the very beginning, because 
the axion dynamics remains unchanged. The property of the 
domain walls, however, may be crucially modified. This is 
the case if the domain wall number N is equal to unity. In this 
case, all the potential minima of Vi along the potential valley 
(02 = 0) are identical, and therefore, one cannot generate any 
bias between them. If domain walls were bounded by strings, 
they would be unstable and collapse quickly after formation. 
In the case of the axion roulette, no strings are formed, and 
so, the domain walls with N = 1 are stable in a cosmological 
time scale [28]. This is nothing but a cosmological disaster, 
and the domain walls must be inflated away. The value of N 
is crucial for the fate of domain walls. 

The domain wall formation implies that the primordial 
quantum fluctuations of the axion result in fluctuations of or¬ 
der unity through the chaotic dynamics during the level cross¬ 
ing. Therefore, at the domain wall formation, the initial en¬ 
ergy densities of domain walls as well as coherent oscillations 
have fluctuations of order unity at superhorizon scales. As 
the Universe expands, the domain wall network is known to 
quickly follow the scaling regime [29-31] Therefore, the en¬ 
ergy density of the domain walls in the scaling regime have 
only fluctuations of order unity at subhorizon scales, and their 
fluctuations at superhorizon scales are suppressed. On the 
other hand, the axion coherent oscillations produced during 
the domain wall formation remain to have density perturba¬ 
tions of order unity at superhorizon scales. Therefore, those 
axions should decay before their energy density exceeds about 
10 “® of the total energy density, since otherwise the non- 
Gaussianity will be too large. If the axion domain walls are 
responsible for the origin of the baryon asymmetry, isocurva¬ 
ture perturbations and their non-Gaussianity are also gener¬ 
ated [12]. In another case, the axion produced by the domain 
wall decays may contribute to the cold dark matter. 

If the axion roulette works, the kinetic energy of the axion 
exceeds the height of the potential barrier, which is impossible 
to realize in the case of a single axion field. Such a large ve¬ 
locity of the axion could be used to enhance the baryon asym¬ 
metry in the spontaneous baryogenesis scenario [9, 10]. We 
leave a detailed analysis of this scenario, especially the esti¬ 
mate of the isocurvature perturbations for future work. 

In summary, the axion roulette works and the axion domain 
walls are formed if (i) the axion starts to oscillate around or 


slightly before the level crossing, and (ii) the initial oscillation 
energy is sufficiently large to climb over the potential barrier 
(cf (7)). The first condition does not require any severe fine- 
tuning; for instance, this is naturally satisfied if the final two 
axion masses are comparable and the decay constants /i and 
/2 are around the GUT scale or higher. The second condi¬ 
tion is met if the alignment mechanism is operative, and it is 
known that such enhancement occurs frequently for several 
axions [32, 33]. Therefore, the domain wall formation may be 
a more common phenomenon in the Axiverse than previously 
thought. 
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